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The trace of the stress-energy tensor of the cosmological fluid, proportional to the Ricci scalar 
curvature in general relativity, is determined on cosmic scales for times ranging from the inflationary 
epoch to the present day in the expanding Universe. The post-inflationary epoch and the thermal 
history of the relativistic fluid, in particular the QCD transition from asymptotic freedom to con- 
finement and the electroweak phase transition, leave significant imprints on the scalar curvature. 
These imprints can be of either sign and are orders of magnitude larger than the values that would 
be obtained by naively extrapolating the pressureless matter of the present epoch back into the 
radiation-dominated epoch. 



I. INTRODUCTION 

The broad history of the Universe, according to the 
Standard Cosmological Model, is a succession of expan- 
sion epochs from the earliest moments of the Big Bang 
to the present day. In reverse chronological order, these 
are: dark-energy domination, during which something 
like a cosmological constant drives the accelerated cos- 
mic expansion; matter-domination, during which cold 
dark matter and baryons shape the large-scale structure 
of the Universe; radiation-domination, during which the 
properties of the cosmic fluid are described by the rela- 
tivistic degrees of freedom of the Standard Model (SM) of 
particle physics and beyond; and inflation, during which 
the potential energy of a scalar degree of freedom dom- 
inates the cosmic energy density, producing an expo- 
nential burst of cosmic expansion. Despite fundamental 
uncertainties — the nature of dark energy and dark mat- 
ter, and the physics beyond the Standard Model remain 
open questions, whilst the search continues for direct 
evidence of the inflationary birth of the Universe — this 
sequence of epochs forms the architecture of a widely- 
accepted, concordant model of the Universe. 

The aim of this article is to give a more quantitatively 
accurate account of cosmic history than one obtains from 
patching together inflationary, radiation-, matter-, and 
dark energy-dominated epochs. It may often suffice to 
describe the cosmic energy density p 
terms of the scale factor a and the equation of state 
w = (—1, 1/3, 0, —1) for inflation, radiation, matter, and 
dark energy. However, it is rare to find a detailed history 
of the trace of the stress-energy tensor of the cosmological 
fluid, or equivalently the Ricci scalar curvature, R = kQ 
where « = 8irG, and 9 = -Tjf = p - ip. 1 The final 



output of our calculations is a determination of R and 
0, as functions of cosmic time or redshift: See Fig. [5J 

To make these calculations, we must properly account 
for the thermal history of the cosmological fluid, the dy- 
namics of the phase transition of quantum chromody- 
namics (QCD), and we must give some account of the 
electroweak phase transition. We do not solve the full 
Boltzmann equations, as necessary for the standard cal- 
culations of the light element abundances in Big Bang 
Nucleosynthesis and the relic dark matter abundance in 
extensions of the Standard Model (e.g. [l|), but instead 
use free field thermodynamics where applicable. We also 
make some assumptions about inflation and reheating in 
the early Universe in order to follow R and all the way 
back to the inflationary epoch. 

It is often supposed that the trace of the stress-energy 
tensor is negligibly small during the radiation-dominated 
epoch. It is also supposed that the dominant contribu- 
tion to the trace comes from baryons and dark matter. 
Both of these suppositions are incorrect: massive parti- 
cle species contribute non-negligibly to the trace when 
the temperature is comparable to their mass (see e.g. 
Refs. 2, 3], wherein the history of O is required to study 
cosmic evolution in a scalar-tensor theory of gravitation) , 
and the effect of phase transitions on the trace of the 
stress-energy tensor is in fact much larger than the small 
contribution from baryons and dark matter. 

Our motivation is to chart a detailed history of as 
comprehensively as we can. Two advances outside of 
cosmology make this possible. First, progress in lattice 
QCD shows that the self-interactions of the quark-gluon 
plasma produce a trace that is very different than what 
one would obtain using free field theory [lj], [l2[ ■ Second, 
the Higgs mass is now known — or at least, the mass of a 
particle whose most straightforward interpretation is the 
Standard Model Higgs [H, E3] • Because we don't know 



1 Conventions: the metric has signature — 
Field Equations are G M " = kT^ v , with R 



+, the Einstein 

^ ~ r ^,A + - and 
The variable is introduced to avoid 



confusing the trace with temperature T. 
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what lies beyond the Standard Model, we offer two sce- 
narios. In the first scenario (despite all arguments about 
naturalness) nothing lies beyond the Standard Model un- 
til one reaches the energy scale of inflation, only some- 
what below the typical scale of Grand Unified Theories. 
In the second scenario, there is a supersymmetric exten- 
sion of the Standard Model plausibly within the reach of 
discovery at the LHC, and then inflationary physics at 
a similar energy scale to before. Our final figure covers 
the entire history of the Universe from inflation to the 
present day for both scenarios. 



II. THE MATTER-DOMINATED EPOCH AND 
AFTERWARDS 

A good (and, by now, entirely conventional) account 
of the expanding universe for z < 100 can be achieved 
by including only non-relativistic baryons, pressureless 
cold dark matter, and dark energy with w = — 1. Of 
course photons and neutrinos are present, but their con- 
tributions to the stress-energy tensor are negligible at 
this stage. As yet, there is no significant evidence for 
time-evolving dark energy (e.g. Ref. [4]); however, the 
inferred, model-independent evolution at z % 1 remains 
the subject of intense investigation as a measure of the 
dark energy equation-of-state. Thus w — —1 for dark 
energy is only a default assumption. 

The trace of the stress-energy tensor, and therefore the 
Ricci scalar curvature, during this epoch is therefore 

R = k<3 = n(p c + p b + (1 - 3w de )pde) (1) 
- 3H 2 ((fi c + n b )(a /a) 3 + An de ) , 

as illustrated in Fig.[TJ The central values of the WMAP- 
derived cosmological parameters — tti>h 2 — 0.02249 for 
baryons, Vt c h 2 = 0.1120 for cold dark matter, and f^e = 
0.727 |5| for dark energy — are used to set pde — QdePo, 
p c = tt c p (a /a) 3 , p b = fl b po(a /a) 3 and p Q = 3Hq/SttG 
where a is the expansion scale factor. The Hubble con- 
stant is H = 100/1 km/s/Mpc with h = 0.738 ± 0.024 
@. For the moment we neglect any non-relativistic neu- 
trino species, whose contribution is small relative to the 
baryons and dark matter. We also ignore all fluctua- 
tions and effects of structure formation: our purpose is 
only to find the Robertson- Walker (RW) background that 
best approximates the expanding Universe at very large 
scales. 



III. FREE FIELD TREATMENT OF THE 
RADIATION-DOMINATED EPOCH 

The free field treatment of the radiation-dominated 
epoch is based on the following standard formulas for 
a collection of non-interacting fermions and bosons in 



thermal equilibrium: 



Piree — 
Pfree — 




where the sum is over all species of particle, nij is the 
mass of the j-th species, gj is its multiplicity, and Sj = +1 
for bosons and —1 for fermions [l|. We have allowed for 
different species to have different temperatures Tj. For 
each collection of mutually thermally equilibrated parti- 
cles, the common temperature must evolve in such a way 
that the collection's contribution to the stress-energy ten- 
sor is conserved. 

It is not necessary to generalize (^J) and ([3]) to include 
chemical potentials, because all the relevant chemical po- 
tentials (for example, for baryon number) are small com- 
pared to the temperature at times when explicit expres- 
sions for non-zero pressure are needed. 

The contribution to the trace O by photons is identi- 
cally zero, whereas a non-zero trace is obtained for any 
free, massive species. In the relativistic limit m <C T the 
leading contribution to the trace is 

w - \^m 2 T 2 boson. 1 ' 

This means that as long as the relativistic fluid includes 
a massive species, then the trace, and therefore the cur- 
vature, is no n- vanishing. The energy density and pres- 
sure decay at slightly different rates as a species cools, 
so that as the temperature drops past about half the 
mass, a bump occurs in the contribution to the trace. 
As an illustrative example, Fig. [5] shows the ratio 8/T 4 
for a single boson or fermion species cooling below its 
rest mass. We expect this process to be repeated for 
each relativistic species of the Standard Model, briefly 
raising the trace above a background level set by any 
abundant, non-relativistic particles for which O = p, or 
else the most massive, relativistic species in equilibrium 
for which 8 oc m 2 T 2 . 

The most recent example of this phenomena is due to 
electrons as the temperature of the cosmic fluid drops 
below the electron mass. Turn ahead to Fig. [6] to see 
this illustrated. Whilst Standard Model neutrinos be- 
come non-relativistic more recently, their contribution to 
the trace is negligible. 2 Consequently, all one needs to 



2 The sum of the neutrino masses is constrained to lie in the range 
0.06 eV < J2 m " < 11 eV at 95 % CL > where the lower bound 
is due to the mass splittings seen in neutrino oscillation exper- 
iments [7j and the upper bound is obtained from cosmological 
observations @, Q. For these masses, neutrinos become non- 
relativistic at a time when the radiation energy density is smaller 
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Figure 1: The Ricci scalar curvature in units of Mpc -2 , or equivalently the trace of the cosmological fluid stress-energy tensor 
in units of GeV 4 , is shown by the thick (blue) curve for 10 -4 < a/ao < 1 or a run of redshift < z < 10 4 . At late times, the 
trace approaches a constant value due to the onset of dark energy domination. For reference, the energy density is shown by 
the thin, long dashed (black) curve. The pressure is shown by the thin, short dashed (red) curve. Note that the pressure is 
negative after a/ao > 0.08, so the absolute value is plotted. In the lower panel, the trace in units of the total energy density is 
shown, 1 — 3p/p. 



model the cosmic trace at temperatures T < 10~ 4 GeV 
or for a/a > 1CP 9 , long before the onset of the matter- 
dominated era, is baryons, cold dark matter, and dark en- 
ergy as given by Eq. |T]). At temperatures T > 10~ 4 GeV 
or a/ao S 10~ 9 , a virtual parade of particle species, each 
cooling past their individual rest mass energies, keeps the 
cosmic trace well above the level of the energy density in 
cold dark matter. 

Proceeding to higher temperatures, it is convenient to 
describe the relativistic gas in terms of an effective num- 



than the contribution due to baryons and cold dark matter. Con- 
sequently, the peak ratio of the trace contributed by neutrinos to 
non-relativistic matter is found to be 0„/© m ~ 0.02 (^2 m^/eV) 
occuring at a redshift z 800 {^2, m v /eV). Hence, this is at most 
a 2% change in the cosmic trace. This is comparable to other 
phenomena which shift the effective number of relativistic neu- 
trino species upwards from 3 to 3.046 which we also neglect. See 
Rcf. 10] for a discussion of these corrections. 



ber of degrees of freedom , g* p for the energy density 
and pressure 

/9 = 5 *30 T4 ' P = 9 * p 90 Ti ' ^ 

The stress-energy trace is 6 = (<?* — g* p )n 2 T 4 /30. For 
the free, massive particles of the Standard Model com- 
posing the cosmological fluid it is straightforward to use 
equilibrium thermodynamics to calculate g*, g* p . As we 
pass to larger temperatures and redshifts, however, the 
main questions are: (1) When can we use a free field 
treatment of thermodynamics? (2) When free field ther- 
modynamics fails, what should we use instead? We will 
consider four cases in the following sections: the QCD 
phase transition, where lattice results are available and 
are significantly different from free field results; the Stan- 
dard Model electroweak transition, where ring-improved 
one-loop calculations indicate that departures from free 
field expectations are not too great; a supersymmetric 
extension of the Standard Model, where we assume that 
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Figure 2: The contribution to the ratio Q/T 4 by a species 
(g = 1) cooling below its rest mass is shown. The solid 
(dashed) line is for fermions (bosons). According to a free 
field treatment, the freeze-out of each species of the Standard 
Model results in a bump to the trace. 



the effect of interactions on the thermodynamics is fairly 
modest; and finally inflation and reheating. 



IV. THE QCD PHASE TRANSITION 

Mesons and baryons melt into constituent quarks and 
gluons at the QCD phase transition. The best source of 
information about thermodynamics in the vicinity of the 
transition is lattice QCD: see for example Refs. (ill. fT3| . 
While the results of these works are qualitatively similar, 
quantitatively they differ enough so that it is best to 
summarize each separately. 

The authors of Ref. [TTJ] provide the following 
parametrized form: 



1 - 



1 



(1 + e (T- Cl )/c 2 )2 



d>2 d± 



(6) 



We will use the values c x = 193.8 MeV, c 2 = 13.6 MeV, 
d 2 = 0.241 GeV 2 , and d 4 = 0.0035 GeV 4 given in [TO], 
which include a 10 MeV shift downward in temperature 
of the lattice data, which slightly reduces the discrepancy 
with [ll| (to be described in the next paragraph). The 
pressure as a function of T can be extracted from the 
general result 



P(T) p{T ) 
2^4 



rp4 
1 



dT 



T 



, e(T') 

y/5 



(7) 



together with some information about p at low temper- 
atures. In order to find p(T) and p(T) at low tempera- 
tures, we used the free field expressions @ and (|3J) with 



the sum running over pions, kaons, and the rj, p, and u 
mesons, together with protons, neutrons, and their anti- 
particles. The free field value of Q/T 4 matches the one 
given in ^ at To = 134 MeV. (It is coincidental that 
this is almost exactly the mass of the neutral pion.) For 
T > Tq we extracted p from ([7]) and then p from (j6|) . 

The authors of [l2[ provide the following parametrized 
form as a fit to their results: 



ho + jo 



tanh(/it + / 2 ) 



1 + 9it + g 2 t 2 



(8) 



where t = T/(200MeV) and h = 0.1396, hi = -0.1800, 
h 2 = 0.0350, f = 2.76, fi = 6.79, g x = -0.47, and 
92 = 1.04. The relation ((7|) can again be used to ex- 
tract p{T). The form © describes the lattice data from 
T = 100 MeV to 1 GeV, and it also is a good fit to the free 
field thermodynamics from T = 20 MeV to T = 100 MeV. 
Thus we can use ([5]) essentially without additional low- 
temperature information — or match to free field thermo- 
dynamics at T = 20 MeV. 

We note that a fairly good fit to the lattice data of 
[12], which covers the range 100 MeV < T < 1 GeV, can 
be achieved with a simpler functional form: 



Q_ 

rp4 



(1 + e (T- Cl )/c 3 )2 



"3/2 
T 3/2 



(9) 



where a = 180 MeV, c 2 = 20 MeV, and d 3/2 = 
0.39GeV 3/2 . The 1/T 3/2 scaling is curious because it 
suggests a leading dimension A ss 5/2 for the opera- 
tor which breaks conformal invariance at higher temper- 
atures. The approximate 1/T 3 / 2 scaling persists over less 
than a factor of 10 on the temperature axis, so we should 
not take this scaling dimension too seriously. No free field 
construction will give a scalar operator with A = 5/2, but 
below 1 GeV we are far from asymptotic freedom. 

The results from [ll| and [l2j both support the conclu- 
sion that the peak in Q/T 4 is narrower and higher than 
one can get from free field expressions: see Fig. [3] Per- 
haps this is not too surprising given that the QCD tran- 
sition is almost a second-order phase transition, blunted 
only weakly into a cross-over, whereas free field expres- 
sions based on a finite particle spectrum are necessarily 
far from describing sharp phase transitions. 

As is evident from Fig. [31 there is still some question 
about just how high and narrow the peak in Q/T 4 is. The 
main difficulty is that there are a number of different lat- 
tice actions, all of which have in principle the same con- 
tinuum limit, but which can in practice give significantly 
different results at finite lattice spacing. Extrapolating 
to the continuum based on finite lattice studies can never 
be fully systematic. The hotQCD results shown in Fig. [3] 
are representative of p4 lattice data with lattice spacing 
a chosen so that N T = 1/ (aT) = 8 in the Euclidean time 
direction. Using instead an asqtad action results in a 
peak roughly 15% lower. The Wuppertal results shown 
in Fig.[3Jare based on a stout-improved staggered fermion 
action with iV r in the range of 8 to 12. (For details on 
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0.2 0.4 0.6 0.8 1.0 
T/GeV 

Figure 3: Comparison of the thermodynamics of lattice QCD 
and of a theory of free quarks and gluons. In all cases, three 
quark species (u, d, and s) have been included. The curve 
labeled "hotQCD" is based on Ref. while the curve la- 
beled "Wuppertal" is from Ref. [l2l ]; see the main text for de- 
tails. The free field theory curve was obtained by setting the 
masses of all the quarks and gluons artificially to 500 MeV. 
This counter-factual assumption is made in order to obtain 
the narrowest possible peak in O/T 4 in the vicinity of the 
QCD transition, illustrating that interactions have a signifi- 
cant effect beyond introducing dynamically generated masses. 

the p4, asqtad, and stout-improved staggered fermion ac- 
tions, see [ll|,ll2| and references therein.) Going forward, 
we will use the p4 results of [ll|; not that we endorse 
them as more accurate, but instead because they are the 
most different from free field results, providing us with 
an upper limit of how sharp the QCD transition might 
plausibly be as a result of interactions. 

In order to obtain a complete account of the visible 
sector degrees of freedom in the vicinity of the QCD tran- 
sition, we combine the energy density and pressure ob- 
tained from ([5]) and ([7]) with free field treatments of all 
the leptons, and also the c and b quark, where all par- 
ticles are constrained to have the same temperature. A 
free field treatment is obviously not perfect (particularly 
for the c quark), but improved approximations would be 
complicated. 

V. THE ELECTROWEAK PHASE TRANSITION 

It is far from obvious that the electroweak transi- 
tion can be approximated by free field thermodynam- 



ics. The masses of all observed Standard Model parti- 
cles owe their existence to a non-zero Higgs expectation 
value, <j) = a = 246 GeV. But this expectation value is 
eventually driven to zero at high temperatures. An elec- 
troweak scale contribution to the cosmological constant 
accompanies this change in the Higgs expectation value. 
Meanwhile, electroweak interactions introduce thermal 
corrections to particle masses. Without accounting prop- 
erly for these thermal corrections and other loop effects, 
the Higgs mass itself, now known to be approximately 
125 GeV in vacuum [13|, would become imaginary 
once the Higgs expectation value falls below the point 
where the tree-level potential is concave up. One of the 
main conclusions of this section, illustrated in Fig. [H is 
that, close to the peak of O/T 4 , free field thermodynam- 
ics based on the vacuum particle spectrum nevertheless 
provides a decent approximation to the ring-improved 
one- loop treatment of [l5[ , which is the simplest account 
of the electroweak transition that avoids obvious inconsis- 
tencies such as imaginary masses. At substantially higher 
temperatures, we will show that the ring-improved one- 
loop treatment predicts a negative value of O/T 4 . 

A more modern understanding of the electroweak tran- 
sition [l6[, based in part on lattice simulations, is that the 
transition is not weakly first order, as predicted by the 
ring-improved one- loop treatment, but is instead a cross- 
over. If anything, we expect the full non-perturbative 
results for the trace O of the stress tensor to be closer 
to the free field results than the ring-improved one-loop 
results are, though it is likely that the ring-improved one- 
loop treatment is still a good guide well above and well 
below the cross-over. 

The treatment of [15j proceeds in three steps: 



1. First one produces thermally improved formulas 
for all the fields using self-energy diagrams. The 
schematic form of these masses is m 2 ((f>,T) = 
m 1ree{ < t ) ) + 9 2 T 2 , where g is a gauge coupling and 
<fi is the Higgs field expectation value. The precise 
forms of all the masses are listed in table HI 



2. Next one assembles an effective potential, correct 
through one-loop order, as follows: 



6 



Vfa T) = V tIee (<j>) + Vg>(<j>, T) + V theim (<t>, T) 
4 KV 



Vtherm(</>, T) 



vac ' 
2\2 



(10) 



m,-(a,T) 2 2"*' 

roc 

y |flr 4 / dxx 2 log(l - s e-v^w^v^) ) 

27r Jo 



-mj (0, T) 4 + 2m, (0, T) 2 m 3 (a, T) 2 - -m, (a, T) 4 



where the sums over j run over the entries in tablelH 
and Sj = ±1 for bosons and fermions, as in ([2]) 

and ([3]). Vvac is the one-loop vacuum correction 
to the effective potential, and its dependence on 
T arises only because, by prescription, one replaces 
the vacuum masses by the thermally improved mass 
formulas. Vtherm(^) T) is the one-loop contribution 
from thermal occupation numbers for all the fields. 

3. At any fixed temperature T, let <j>(T) be the value of 
the Higgs field which minimizes V((f>, T). Then the 
pressure is p = —V(</>(T),T). The energy density 
can be extracted from the first law formula p = 
-p + Tdp/dT. 




T (GeV) 



Figure 4: Comparison of the ring-improved one-loop treat- 
ment of electroweak thermodynamics [HI] with free field ex- 
pectations, shown by the solid black and dashed red lines, 
respectively. Both curves are based on the minimal Standard 
Model with the Higgs of mass m H = 125 GeV [11 03]. The 
sharp feature at 160 GeV, corresponding to the temperature 
at which the symmetry of the Higgs potential is broken, is a 
shortcoming of the ring-improved one-loop treatment. 

In Fig. 21 we compare the results of the ring-improved 
one-loop approach to the case of free field thermodynam- 
ics based on the masses and multiplicities of the Standard 
Model particle spectrum at zero temperature. The kink 
in O/T 4 at 160 GeV shows the location of the first order 



transition predicted by the ring-improved one-loop ap- 
proach. This kink is actually a very slight discontinuity 
in O/T 4 as a function of T. The kink is smoothed out by 
the non-perturbative effects that make the entire transi- 
tion a crossover. It is clear from Fig. |4] that electroweak 
interactions and symmetry restoration have essentially 
no effect on the peak value of 0/T 4 , and that they are 
mainly significant on the high-temperature side of this 
peak. In the next two paragraphs we continue our inves- 
tigation of the high-temperature behavior analytically. 

At temperatures high above the electroweak transition, 
we expect that the leading contribution to the trace is 
due to the most massive relativistic species in the ther- 
mal bath. In the electroweak symmetric phase, however, 
the masses of most particles vanish. Although the longi- 
tudinal gauge bosons jz, Wl 7 Zl have large, thermally- 
improved masses, these temperature-dependent correc- 
tions drop out when we evaluate p—Zp = —Ap + Tdp/dT. 
This is because their thermally-improved masses, in the 
high-temperature, symmetry-restored phase, are dimen- 
sionless constants times the temperature, resulting in a 
contribution to p that scales precisely as T 4 . The only 
particle in the spectrum whose mass in the symmetry- 
restored phase is not linear in the temperature is the 
Higgs. So its thermal excitations dominate O at high 
temperatures. An important sub-dominant contribution 
is the classical Higgs potential. Altogether, at T 3> a, 

9 = h 2 T 2 + h a + 0(T- 2 ) , (11) 

where 

h = \^J 2 m 2 H + (one loop) (12) 

h% = — ( — —J / dx — ^ = ^ = j^=^ 

V 2tt J J ^ x 2 +(?2 ( e V* 2 +« 2 _ i) 

q 2 = (m 2 H + 2m 2 vac + 2m 2 Vvac + m 2 Zvsic )/(2<T) 2 . 

The contribution to ho shown explicitly is from the tree- 
level Higgs potential. The corrections, indicated by 
(one loop), are numerically small given the measured 
value of the Higgs mass. The h^T 2 term is from the 
thermal excitations of the Higgs. Numerically, — hi ~ 
m 2 ? /45 ~ (19 GeV) 2 , valid for T > a. In the absence 
of contributions beyond the Standard Model, the trace 
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Particle 


Mass squared 
(m 2 ) 


Multiplicity 
(3) 


W 


771 VK, vac ^7 


4 


W L 


2 / 2 , 22 T 2 \ 
m W,vac T — J 


2 


7 





1 


v 
Z 


m Z,vac^7 


z 


7i 


\{9i+9i)C{T) 


1 




Z L 


\(gi + gi)C{T) 


1 


C(T) = \<t> 2 + fT 2 


H 


Ao(3^ 2 - ( 7 2 ) + iAoT 2 
+ 3 72-(2m 2 vac + 2mfv iVa[: + m|_ vac ) 


1 


t 


m t,vac ~7 


12 


other 
fermions 





78 


gluons 
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Table I: Thermally corrected masses of Standard Model par- 
ticles as a function of the Higgs expectation value <j> and the 
temperature T, following [Hjj |. A few entries require explana- 
tion. At T = 0, 7l becomes a physical, massless photon, but 
at T 7^ it mixes with the Higgs scalars along with the lon- 
gitudinal Wl and Zl ■ The H entry corresponds to quanta of 
the magnitude <f> of the Higgs doublet. We ignore the masses 
of all the fermions except the top quark, and we also ignore 
strong interactions. 



Quantity 


Value 


Quantity 


Value 


a 

9i 


246 GeV 
0.357 


Ao 

92 


2a' 2 

0.652 



Table II: Standard Model parameters. 



is negative at high temperatures, starting at about 
the temperature T* = y/— h a /h 2 ~ 1200 GeV where 
the negative h-fT" 1 term in (|TT|) overtakes the positive 
ho term. The most negative value attained by 0/T 4 

1600 GeV at which point 
—2 x 10 -4 . Because this 



occurs at T = d— 2ho/h2 - 
9* - g*. P = -I5hl/(2ir 2 h ) 
negative trace is a result found using the ring-improved 
one-loop approximation, it is not to be regarded as iron- 
clad. We leave for future work the interesting question 




10 10 
T/GeV 

Figure 5: The degrees of freedom of the Standard Model and 
a supersymmetric extension are illustrated as functions of 
equilibrium temperature. In the top panel, the degrees of 
freedom of the energy density, g, (heavy curves), and pres- 
sure g* p (light curves), are shown for the Standard Model 
(solid curves) and a representative supersymmetric extension 
(dashed curves) . In the lower panel, the difference g t — g tp is 
shown. Note that adding a single dark matter species to the 
Standard Model curves does not produce a discernible change 
in these figures. 



of whether it persists beyond this approximation. 

We may reasonably expect the LHC or direct dark 
matter detection experiments to discover new particle 
species at higher masses. Consequently, we supplement 
the Standard Model particle spectrum with a dark mat- 
ter (DM) particle. Until it freezes out, a single, relativis- 



tic, bosonic species contributes <d 



DM 



m 



DM 



T 2 /12 to 



the trace. If a bosonic dark matter particle has a mass 
m DM < y/l2h 2 ~ 70 GeV or m DM > 20TeV then there 
will be a window in which the collective stress-energy 
trace is negative due to the Standard Model Higgs. Such 
a possibility is reflected in Fig. [6l where we include a 
100 TeV dark matter particle which freezes out after it 
becomes non-relativistic, in our extrapolation of O back 
to times before and temperatures above the Standard 
Model electroweak transition. 

In hindsight, it is perhaps not too surprising that a 
weakly coupled theory like the Standard Model manages 
to have thermodynamic properties not too far from free 
field theory. We therefore feel justified to include here 
a similar free field account of the constrained Minimal 
Supersymmetric Standard Model (cMSSM) , in which we 
use DarkSUSY [17] with input parameters correspond- 
ing to the best fit to the recent LHC data (Table 4 
of Ref. [3: m = 389.51 GeV, m 1/2 = 853.03 GeV, 
A = -2664.79 GeV, tan/3 = 14.5, and sgn(^) = +1) 
to generate a zero-temperature mass spectrum, which we 
then plug into © and ^ without regard to the mod- 
ifications that will arise from interactions and symme- 



8 



try restoration. Though these modifications are proba- 
bly significant on the high-temperature side of the peak 
in 0/T 4 , the main uncertainty lies not in their precise 
determination, but our ignorance of what the superpart- 
ner spectrum really looks like (assuming it is there at 
all). As seen in Figs. [5] and [6j an anomalous trace arises 
as the spectrum of sparticles with masses in the range 
~ 0.3 — 2 TeV become non-relativistic. More specifi- 
cally: the best-fit model we used as a benchmark has a 
neutralino LSP at 363 GeV, additional color singlet spar- 
ticles below 1 TeV, and colored sparticles above 1.4 TeV 
except for a stop at approximately 817 GeV. 

Proceeding to ever higher energies requires a theory 
to guide our calculation of the trace. In the approach 
to GUT-scale energies, we may pass through a desert, or 
there may be new particle species littering the landscape. 
We should also expect couplings to run, although this is a 
small effect in simple extensions of the Standard Model. 

Pausing to reflect on Fig. [6j we may assess the degree 
of error commonly made in the literature by assuming 
the curvature to be negligible, or else solely determined 
by the non-relativistic matter content extrapolated back 
from the matter-dominated era. Instead, we see that at 
various epochs, when the number of degrees of freedom 
changes sharply, the trace rises to an appreciable fraction 
of the total energy density. 



VI. INFLATION AND REHEATING 

The history of curvature at earlier times requires a the- 
ory of the early Universe, which we take to be the infla- 
tionary dynamics of a scalar field. Whereas other theories 
of the primordial Universe will yield different predictions 
for the trace and curvature, the merit of inflation is the 
consistency with observations and ubiquity of examples. 
The generic course of inflation is exponential expansion 
proceeding with pi ~ —pi caused by the slow roll of the 
inflaton i\) down its potential until some late stage when it 
breaks from slow roll and plummets down the potential 
towards the vacuum. In broad terms, the inflaton un- 
dergoes damped evolution as it decays into a turbulent 
mixture of relativistic particles and massless radiation. 
It is only after some transient period that these decay 
products equilibrate and can be appropriately described 
as a thermal bath with temperature T. The description 
of the cosmological fluid, trace and curvature, can then 
merge with our description of the high temperature limit 
beyond the Standard Model. 

A detailed analysis of the post-inflationary, preheat- 
ing epoch including the approach to equilibrium requires 
a numerical simulation of the nonlinear, quantum fields 
{e.g. Refs. [l9l - l2ll ]). Nonetheless, we can capture the 
basic behavior with a semi-analytic model. To be defi- 
nite, we consider an inflaton with potential V = ^m 2 ijj 2 , 
which decays into relativistic particles and massless ra- 
diation via a damping term. Following Refs. P, we 



evolve the system 

i) + (3H + r)ip + V' = 0, (13) 
Pr + AHpr-Tij 2 = 0, 

with pi = \i? + V, pi = \ip 2 - V. The field is started 
in a slow-roll state, with m ~ 10 GeV chosen to yield 
satisfactory inflation fFi 2~> . As slow roll ends, the in- 
flaton begins to oscillate about the potential minimum, 
such that pi — ipi swings between +4pi and — 2p.;; despite 
the fact that the time-averaged inflaton pressure is small 
compared to the energy density, the trace and curvature 
become negative. This pressureless phase is shortlived, as 
an assumed inflaton decay width T ~ 10 10 GeV causes 
the field to decay exponentially with cosmic time; the 
trace O; = pi — Zpi plummets until it meets the anoma- 
lous trace p r — 3p r cx m 2 T 2 due to any massive species 
among the inflaton decay products. This scenario, be- 
ginning from the end stages of inflation, is illustrated in 

Fig.m 

The variety of inflationary models, and the challenges 
of accurately simulating the nonlinear dynamics of re- 
heating, translates into a wide degree of uncertainty in 
the behavior of the trace and curvature in the approach 
to equilibration. In the above model we have assumed 
an idealized description of the radiation produced by the 
decay of the inflaton. If the inflaton decay products do 
not achieve the relativistic equation of state w = 1/3 un- 
til after a brief equilibration period, then the trace may 
be much larger. If one of the decay products includes 
a kinetic-energy dominated scalar field with equation of 
state w — +1 then the trace contains a residual, nega- 
tive term that decays as cx a -6 . This possibility is also 
illustrated in Fig. [7] 

The gap between the trace and the energy density is 
so wide in the interval between the end of inflation and 
the electroweak epoch, assuming no new physics, that 
we may question whether there are any other sources of 
curvature or trace. For example, long-wavelength fluc- 
tuations in the energy density and pressure might lead 
to an rms contribution. If the pressure and energy den- 
sity perturbations imparted by inflation unto the cosmo- 
logical fluid are characterized by a propagation speed v 
where v 2 = Sp/Sp, then our question becomes: is v 2 as 
closely tuned to 1/3 as is w, the background equation of 
state? In the case of perfectly adiabatic perturbations, 
with v 2 = w, then the answer is yes. At some level, how- 
ever, we expect a departure from adiabaticity. Recall 
that the linearized perturbation to the Ricci scalar in a 
coordinate system 

ds 2 = -dt 2 + a 2 {t)e 2 ^ t '^8 ij dx i dx j (14) 

is SR = 25R^ - 6(C + 4H"C), and SR^ = 2V 2 C/a 2 , 
where £ is a familiar variable from studies of inflation. 
(See Ref. 26].) For perturbations of wavelength much 
greater than the radiation-era Hubble radius, specifically 
those perturbations responsible for the observed large- 
scale structure, then ( is a constant in time and SR^ 
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Figure 6: The Ricci scalar curvature in units of m 2 , or equivalently the trace of the cosmological fluid stress-energy tensor 
in units of GeV , as in Fig. 1, for 10" 18 < a/a < 10" 4 . The energy density and the pressure are also shown as long dashed 
(black) and short dashed (red) lines. Except as noted below, = p — 3p is positive. 

For reference, a continuation of the energy density in pressureless matter is shown by the dotted black line; naive treatments 
that equate the trace with the matter density at early times are easily off by ~ 5 — 10 orders of magnitude. The solid (blue) 
curve shows for the Standard Model plus a dark matter particle of mass tudm — 100 TeV. The thin dashed (blue) curve shows 
the case of the SM through the electroweak transition as described in the text, with becoming negative at a/ao < 10 -16 . 
The long dashed (blue) curve near the left edge of the plot shows the case of the cMSSM. In the lower panel, the trace in 
units of the critical density is shown; the spikes at a/ao = 10~ 16 , 10~ 15 , 4 x 10~ 13 , 10~ 9 correspond to particles of the cMSSM 
becoming non-relativistic, the electroweak transition, the QCD transition, and the departure of electrons from equilibrium as 
they become non-relativistic, respectively. 



becomes negligibly small so that 5R -C R- In this case, 
the homogeneous quantities we have calculated provide 
the leading contribution. 

To conclude this brief history of curvature, we show all 
three epochs in a single display, in Fig. |SJ 



VII. DISCUSSION 

The history of the curvature and trace are of interest, 
beyond their intrinsic value as a record of the significant 
events in the large-scale evolution of our Universe. For 



example, an accounting of the thermal history of the uni- 
verse is necessary to accurately account for the spectrum 
of inflation-produced gravitational waves [27[. However, 
the most dramatic effect may arise when considering new 
gravitational phenomena beyond general relativity. 

Gravitational theories that introduce new degrees of 
freedom may require an accurate model of the trace 
of the stress-energy tensor. The canonical example is 
the nonminimally-coupled scalar field, for which the La- 
grangian is 

C = i (1 ~ 2 ^) R ~ \^) 2 - v i<P) + £sm, (15) 
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Figure 7: The Ricci scalar curvature in units of fm -2 , or equivalently the trace of the cosmological fluid stress-energy tensor in 
units of GeV 4 , as in Fig. [TJ for 10~ 32 < a/ao < 1CF 18 . Except as noted below, = p — 3p is positive. The energy density and 
the absolute value of the pressure are also shown as long dashed (black) and short dashed (red) lines. 

Inflation continues at earlier times than shown in the plot, with all curves continuing as horizontal lines. The solid (blue) line 
shows the extrapolation of the SM with a DM species with mass 100 TeV back to the post-inflationary epoch. The long dashed 
(blue) line shows the extrapolation of the cMSSM where the LSP is at 363 GeV. In the cMSSM scenario, a putative decay 
product of the inflaton having equation-of-state +1 is considered, for which the contribution to the trace is negative. In such 
a case, the curvature and trace may become negative, as shown in the range 3 x 10~ 27 < a/ao < 3 x 10~ 23 by the thin dotted 
(blue) line. also becomes negative during oscillations at the end of inflation. In the lower panel, the trace in units of the 
critical density is shown. After the inflaton makes several oscillations across the bottom of its potential, and until the Universe 
is fully radiation dominated, we plot the time average of the trace for clarity of display. 



and where Csm includes the Standard Model. The ip 
equation of motion is dtp = V + 2£Rip so the curvature 
introduces a novel, time- and space-dependent effective 
mass such that when £R < then the field is susceptible 
to a tachyonic instability (e.g. [lllH^]). The curvature 
is set by the dynamics of the field, which are in turn 
sourced by the trace of the stress-energy of any other 
forms of matter. Specifically, 

1 - 2k 2 ^ 2 (1 - 120 
In the case of a strongly coupled field with |£| ^> 1, for 



which \£k<p\ <C 1 so that the ip field does not dominate 
the stress-energy, then we have approximately R ~ kO. 
In the case of a conformally-coupled field, £ = 1/12 and 
V = \<p 4 , then R = kQ, which means that the trace cal- 
culated in this paper is directly applicable to the evolu- 
tion of <p: every feature in our O curve results in a feature 
in ip. Similarly, nonlocal quantum gravitational effects 
that depend upon quantities such as D _1 i? or A" 1 ./? 2 , 
where A p is a fourth-order conformally-invariant oper- 
ator [lO, [Ml, are sensitive to the history of curvature. 
Consequently, the results presented in this paper for the 
evolution of the trace may prove useful in detailed stud- 



ies of new gravitational physics, such as gravitationally- 
coupled inflation |32j or dark-energy motivated modifi- 
cations of gravity (e.g. Ref. [33| . or [34| and references 
therein). 
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Figure 8: A cosmic history of the Ricci scalar curvature in units of m -2 , or equivalently the trace of the cosmological fluid stress- 
energy tensor in units of GeV 4 . The run of time scales is from the end stages of inflation to the present day, 10 -32 < a/ao < 1. 
The energy density and the absolute value of the pressure are shown by the thin long dashed (black) and thin short dashed 
(red) curves. Also shown is an extrapolation of the energy density in pressureless matter into the radiation era, for reference, 
which terminates near the possible freeze-out of a dark matter particle of mass mr>M = 100 TeV. This plot is a composite of 
Figs. [1] [6] and [7] and the explanatory comments of previous captions apply to this figure as well. 

In the lower panel, the trace in units of the critical density is shown. The inset box is magnified by a factor of three to show 
spikes in curvature, due to the detailed thermal history of the SM and a supersymmetric extension, the electroweak transition, 
the QCD transition, and the departure from equilibrium of electrons. Simple treatments that attribute the curvature due to 
the contribution to the trace by the energy density in pressureless matter at early times underestimate the amplitude of the 
trace by many orders of magnitude, and miss out on episodes in which the curvature and trace are negative. 
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